
FIGURE 1. Sketch of colliding flux tubes in cross section. The solid lines are the tube 
boundaries. The radius of the isolated tube is R, the length of the interaction region is L, 
the y axis is chosen parallel to the contact surface, and e is the fractional deformation. 
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FIGURE 2. (a) Sketch of the magnetic flux surfaces of the colliding tubes, projected onto 
the axial cross section. Reconnection can take place within a thin layer along the contact 
surface, (b) Sketch of the Parker-Sweet reconnection geometry, which we assert ensues once 
the configuration shown in Figure (2a) is established by the collision. The y axis is again 
the contact surface and V{ n corresponds to the incoming tubes. In slow collisions omly a 
slight pileup of magnetic flux is sufficient to permit Sweet-Parker reconnection. 
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FIGURE 3. Plot of the function jfo /4 J(y ) as a function of y over its permitted range. As 
eqn. (3.12) shows, this function is proportional to the amount of magnetic flux reconnected 
in a collision. 
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ABSTRACT 

We develop an analytical theory of reconnection between colliding, twisted magnetic flux 
tubes. Our analysis is restricted to direct collisions between parallel tubes and is based 
on the collision dynamics worked out by Bogdan (1984). We show that there is a range 
of collision velocities for which neutral point reconnection of the Parker- Sweet type can 
occur, and a smaller range for which reconnection leads to coalescence. Mean velocities 
within the solar convection zone are probably significantly greater than the upper limit 
for coalescence. This suggests that the majority of flux tube collisions do not result in 
merging, unless the frictional coupling of the tubes to the background flow is extremely 
strong. 

1. INTRODUCTION 

Observation, theory, and numerical experiment all suggest that magnetic fields become 
concentrated into sheets and ropes when embedded in gas pressure dominated, turbulent, 
highly conducting fluids such as stellar convection zones (Galloway & Proctor 1983, Schmitt 
& Rosner 1983, Hughes and Proctor 1988) and possibly accretion disks (Stella & Rosner 
1984, Sakimoto & Coroniti 1989, Schramkowski & Achterberg 1993). These flux tubes 
move semi- independently of the surrounding fluid: accelerated by buoyancy and magnetic 
tension forces but coupled to the ambient medium by frictional drag (e.g. Parker 1982, 
Moreno- Insertis 1986, Choudhuri & Gilman 1987, Chou & Fisher 1989, Fan, Fisher, & 
DeLuca 1993, DeLuca, Fisher, & Patten 1993). 

Under these conditions, magnetic flux tubes are expected to collide with one another. 
Collisions might result in coalescence, and possible reconnection of the magnetic field, or 
the tubes might retain their identities and separate. It is clear that the rules which govern 
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collisional interaction, and in particular the reconnection process, play an important role 
in the evolution of a flux tube dominated magnetic field. 

Bogdan (1984; hereafter B84 ) considered collisions between magnetic flux tubes and 
developed a model for their merging. He took the flux tubes prior to the collision to 
be infinite cylinders with both axial and azimuthal components of magnetic field, and he 
assumed that during the collision the tubes would flatten along a common surface. If the 
tubes are twisted in the same sense, their azimuthal components in the interaction region 
are antiparallel, and might be expected to reconnect. B84 treated the reconnection as a 
ID problem in magnetic diffusion across the common boundary (with the thickness of the 
boundary taken as a free parameter), and he assumed that the tubes coalesce if the energy 
of the reconnected magnetic field exceeds the relative kinetic energy of the collision. The 
amount of reconnected flux is proportional to the magnetic diffusivity rj which is typically 
small in situations of interest, so coalescence occurs only for very slow collisions. 

The purpose of this paper is to reconsider magnetic reconnection within the framework 
of B84 's collision model. We argue that reconnection is more likely to proceed through a 
magnetic x— point and to be of the Parker-Sweet type. This determines the width of the 
diffusion layer and leads to a scaling in which the reconnection rate is proportional to 772 . 
This results in a model in which flux tubes can coalesce even when moving much faster 
than the critical velocity for coalescence estimated by B84 . 

In Section 2 of this paper we review the dynamics of flux tube collisions. In Section 
3 we describe the reconnection phase and in Section 4 we estimate the critical velocity 
below which two flux tubes can merge. Section 5 is a critical review of the assumptions, 
application of the theory to the solar convection zone, and discussion. 

2. COLLISIONS BETWEEN MAGNETIC FLUX TUBES 

We consider head-on collisions between two identical, parallel magnetic flux tubes. 
Following B84 , we assume that before contact the tubes are circular cylinders with radius 
R and that each tube is in magnetostatic equilibrium, with a helical magnetic field zB z (r)-\- 
(pB^r) and gas pressure -P(r), where r, z, and <f> are cylindrical polar coordinates for each 
tube. We also assume that the relative velocity is slow enough that the flows inside and 
outside the tubes are incompressible. The medium surrounding the tubes is unmagnetized 
and has uniform density p e . 

The equation of motion of the tubes prior to contact was solved more than a century 
ago by Hicks (1879), and the solution is reviewed by B84 . The mass per unit length M. of 
each tube is increased to an effective mass M e ft due to the entrainment of ambient fluid: 
M e f£ — M- is of order 7rp e R 2 . In stellar interiors the ratio of magnetic to gas pressure inside 
the tubes is probably rather low, so p e is probably nearly the same as the density p in the 
tubes, and Jv ^ t should be of order unity. 

If the tubes approach each other with speed \voo in the center of mass frame, the 
actual relative velocity v c at impact tends to be reduced from by drag effects, but can 
also be increased by center of mass motion perpendicular to the direction of approach, 
and decreased by center of mass motion parallel to the direction of approach. All of these 
effects are small for M e ft ~ M.-, and we will parameterize the collisions by the relative 
velocity at impact, v c , and refer the reader to eqn. (18) of B84 for the exact relationship 



2 



between v c , Uoo, and the velocity of the center of mass. 

When the tubes collide, they become flattened along a common face. We measure 
the deformation of each tube by the displacement of a boundary point on the shortest line 
connecting the z— axes of the tubes, and denote this quantity by Re(t). Then (see Figure 
1), the length L of the common face is 

L = (8e)%R. (2.1) 

The deformation of the tubes is maintained by an overpressure AP which decelerates 
the tubes and eventually forces them apart. B84 showed that for small deformations 
AP oc ePj(P), and wrote AP as a power series in the distance | y | from the midpoint of 
the common boundary (see Figure 1) 

The coefficients p n were left undetermined by B84 , so we treat them as unknowns. 

Since the overpressure acts along the length L of the common boundary, eqns. (2.1) 
and (2.2) show that the restoring force per unit length of tube is proportional to ez. The 
equation of motion for e is 

The problem up to this point has an intriguing resemblance to the problem of encounters 
between finite, isotropic elastic bodies (Landau and Lifshitz 1970), which has been solved 
completely and for which it is also found that the repulsive force between the bodies is 
proportional to the | power of the deformation parameter. We therefore thought that we 
might be able to find the coefficients p n by developing an analogy with elasticity theory. 
However, the force between infinite elastic solid cylinders is found to be proportional to 
the first power of the deformation parameter, and we were unable to completely cast the 
MHD equations in the form of elasticity theory (see Appendix), so the analogy proved 
fruitless and the p n remain undetermined. 
Equation (2.3) has the first integral 

l ( de 2 ^Bl(R) ^ Pn f j£_ 

2W 57rM eff ^(2n + l) 8P 2 ' K ' ) 

where the constant of integration follows from the definitions of e and v c . We now convert 
eqn. (2.4) to a form that is more convenient for the analysis later in the paper. We 
introduce the Alfven velocity va defined with respect to the azimuthal magnetic field and 
mass density at the edge of the tube 
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define a parameter // which measures the ratio of M e ft to the mass per unit length the 
tube would have if its density throughout were equal to the density at the edge 



_ M e « 



and let the effective Alfven Mach number for the collision be 

,2 



A 



We expect that // is of order unity and Mac is small for the high /3, slowly moving tubes 
in the solar convection zone. Finally, we subsume the unknown coefficients p n into a 
dimensionless parameter 



a^-Y-^-. (2.8) 



Equation (2.4) then becomes 



<f)'-^(*3 c -«*>. (^) 



from which we see that the maximum deformation parameter e max is given by 



emax = (^)^. (2.10) 

a 

Since a is expected to be of order unity, if we restrict ourselves to small deformations then 
Mac must be small. The deceleration of the tubes is gentle until e is close to e max ; v has 
dropped to 0.5u c at e ~ 0.89e max , and to 0.1u c at e ~ 0.996e max . 

If the flux tubes are twisted in the same sense, the azimuthal magnetic field compo- 
nents at the deformed interface oppose each other. This suggests the possibility of magnetic 
reconnection. We estimate the conditions for and amount of reconnection which occurs in 
the following section. 



3. MAGNETIC RECONNECTION 

3.1 Solutions with Flux Buildup 

Reconnection between colliding magnetic flux tubes is treated as a diffusion problem 
by B84 ; the azimuthal magnetic field components in the two tubes are parallel and an- 
tiparallel to the contact surface and diffuse in one dimension across a layer which varies in 
width over time. The reconnection rate derived by B84 is linear in the plasma resistivity rj 
and involves an unknown lengthscale: the width of the layer across which the fields diffuse. 

We pursue an alternative model of reconnection. It appears likely, considering the 
time history of the collision, that a magnetic x— point forms at the initial point of contact 
between the tubes. The situation at later times then resembles that sketched in Figure 
2a. This configuration is strongly reminiscent of Parker-Sweet reconnection (Parker 1957, 
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Sweet 1958a, b), and the reconnection is expected to be much faster than 1 — D diffusion. 
In this section we find the conditions under which Parker-Sweet reconnection can occur. 

We first briefly review the reconnection model (see Parker 1979a for a more complete 
treatment). Consider the magnetic field shown in Figure 2b. The characteristic length 
of the system is L. Resistivity is important only near the x— point, in a layer of width /. 
Fluid is ejected along the y axis at the Alfven speed va = v —- The flow is assumed to 

(47rp) 2 

be incompressible, and therefore the speed Vi n at which material flows toward the y axis 
is related to va by 

Lv in = lv A - (3.1) 

The electric field E z is determined in the diffusion region (near the x— point) by resistivity, 
and in the ideal region (everywhere else) by inductive effects. Since E z must be spatially 
constant in a steady state, we have 



rjB_y 
I 



v in B y , (3.2) 



where the magnetic diffusivity rj = for electrical conductivity a and we have approxi- 
mated the current density by 

J ~ . 

4vr/ 

Equations (3.1) and (3.2) together imply that 

fv A r]\ \ _ 

Vm= \~~[~) = v ps ■ (3-3) 

Defining the magnetic Reynolds number R m = — we have vps = vaRtu 2 and / = 
_ i 

LR m 2 = lps- 

In applying this reconnection model to colliding flux tubes, we identify B y with B^R) 
and p with its value near the edge of the tube. According to the results of Section 2, the 
dynamics of the collision between the flux tubes determines the relationship between the 
collision velocity, Rj^ = e and the length of the contact region, L = (8e) 2 R. Therefore, it 
would appear that eqn. (3.3) can be satisfied only under special circumstances. 

Suppose the incoming velocity Re < vps- In that case, advection of the field by 
the incoming flow can be balanced by diffusion in a layer of width / % ^PS^f~- O n the 
other hand, if Re > Dps, the fieldlines are being brought in too rapidly to reconnect and 
magnetic flux piles up in the layer. Flux pileup has two effects: va is increased, thereby 
increasing vps (see eqn. (3.3)), and the incoming fluid is decelerated by the magnetic 
pressure gradient. The fluid velocity and vps adjust until they balance, and thus we reach 
a state in which Parker-Sweet reconnection can take place. Note that flux pileup is not 
inconsistent with the assumption of incompressible flow, because the fluid can adjust its 
density by motion parallel to the common boundary without perturbing the magnetic field. 
We can expect this to occur as long as the tubes are gas pressure dominated or weakly 
twisted. 
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We quantify these ideas and derive the flux pileup behavior by using conservation of 
momentum to relate the magnetic field B ( / > (R) at the outer edge of the deceleration layer 
to the magnetic field B^i at the inner edge of the deceleration region, which is also the 
edge of the reconnection region: 

/E>-\2 , 2 , B ll /o A \ 

p( Re ) + — £ = PV in + , ( 3 - 4 ) 

Dividing eqn (3.4) by p, substituting for Vi n using eqn. (3.3) with va set equal to vau 
the Alfven speed at the inner edge of the deceleration region, and introducing the Ohmic 
diffusion velocity 

v D = - = — , (3.5) 

Li l\ m 

we derive a quadratic equation for vai with solution 

vai = (v 2 A + 2(Re) 2 +v 2 D )% -v D . (3.6) 

If the collisions are slow and R m is large, the magnetic field is only slightly increased across 
the deceleration layer. To see this, recall that Re < ^v c and therefore is of order Mac^a- 
Using this result we see from eqn. (3.6) that 

^(l + ^MlcHR-^-R- 1 . (3.7) 

V A 

Equation (3.7) holds only if ^(R) ^ ^' 

In summary, consider the history of a flux tube collision from first contact to maximum 
deformation. According to eqn. (2.9), the speed of each tube decreases monotonically from 

to zero. Initially, vps is large because the length of the common boundary, L, is small 
(see eqn. [3.3]), and fieldlines reconnect only through diffusion. But vps drops as the 
contact surface lengthens, and so a time may come when Re > vps- This marks the onset 
of the Parker-Sweet reconnection phase. This phase ends when Re drops below vps once 
again, due to the deceleration of the flux tubes. At this point the system is not driven hard 
enough for Parker-Sweet reconnection to continue, and any residual reconnection occurs 
only through Ohmic diffusion. 

If the collision is sufficiently weak, Re never exceeds ups, and the Parker-Sweet re- 
connection phase is never established. We can establish a criterion for the existence of 
a reconnection phase by substituting eqns. (2.1) and (2.9) into eqn. (3.3) to yield the 
condition that the tube velocity equals vps- 

eHe max -e*)= * , (3.8) 
2 2 aR mg 

where R m g = iJ m (f ) 1S the global magnetic Reynolds number. There is a reconnection 
phase only when the equality (3.8) can be fulfilled, which is seen to be when 

/ 6 s a i 

€max > (— ) 3 =e c - (3.9) 

2 2 5aR mg 
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Or, eliminating e max in favor of Mac using eqn.(2.10) and introducing a global Ohmic 
drift velocity 



_ V A 
VDg = — , 

we find that Parker- Sweet reconnection occurs for 

2»32 a i -L a i _l 

> -^(-)-^A 2 ^ 2 5 ~ 1.65(-)^v^v^ =v crit . (3.10) 
512 /i a /i 

Equations (3.10) shows that an approximate criterion for the existence of a Parker-Sweet 
reconnection phase is that the relative collision velocity exceeds the geometric mean of 
the global Ohmic diffusion speed and the azimuthal component of the Alfven speed. This 
result is almost independent of the unknown parameter a introduced in eqn. (2.8) and 
the parameter fx (introduced in eqn. (2.6)) which is determined by the various quantities 
characterizing the flux tubes and background medium. We expect V£>g C «i in a highly 
conducting medium, so there is a wide range over which the collisions are rapid enough to 
force reconnection but slow enough to be described by the small Mac approximation. 

3.2. The Amount of Reconnected Flux 

We now calculate the amount of magnetic flux that is reconnected during the Parker- 
Sweet phase. Let the two positive, real roots of eqn. (3.8), which delimit the interval 
during which reconnection occurs, be and ej, and let the corresponding times be ti and 
tf. The reconnected flux per unit length in z is 

$ R = f dtB^vps, (3.11) 
Jti 

where B^\ and vps are both functions of time and are determined by eqns. (3.3) and (3.4). 
Equation (3.7) shows that we can set B^\ ~ B^R). Equation (3.11) can be rewritten as 
an integral over e using eqns (2.1), (2.8), (2.9), and (3.3). We find 

®R M AC 1/4. ^ 

- y i{yo), ( 3 - 12 ) 



where y 



RB^ a 2 / 5 



" Vs dv 

i(yo) = I Bsl , (3-13) 



and the parameter yo is defined as 



yo 



vi y 4 (i-y 2 ) 2 



li 2 a 2 / 5 

° Jx mg lvl AC 



and m and yj are the positive real roots of the equation 



y*-y 3 -yo 2 =o, (3.14) 
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which follows from dividing eqn. (3.8) by e^ nax . It follows from eqn. (3.9) that yo cannot 
exceed y ,max = tj. 

6 5 

The integral /(yo) and the limits of integration m and yj cannot be found analytically, 
but they can be found approximately in certain cases, and this is useful in developing a 
scaling for the reconnected flux. Near threshold, when yo is slightly less than j/o,max, we 
find 

9 c 1 

4 • 6 io o 2 

i(yo) ~ — -3 — , (3-15) 

5 2 

where 8 = yo,max — yo- For tube velocities much greater than the threshold given by eqn. 
(3.10), yo approaches zero and m % 0, yj % 1. In this case 

/(yo = 0) = I max = 2 7iffi ~ 1 - 822 - ( 3 - 16 ) 

5r (5) 

We have computed numerically as well. Figure (3) is a plot of y ( J / ' 4 /(yo) versus yo. 

As eqn. (3.12) shows, this quantity is proportional to the reconnection rate. Note that 

there is a maximum value of when Mac is near its threshold value for reconnection 

increases with increasing Mach number, primarily because the interval during which 

reconnection can take place is increasing. As Mac increases further, reconnection can 

take place at almost all times in the collision, but the collision itself takes less time (the 

_ i 

collision time scales as M A q), so begins to decrease again. 

4. THE CRITICAL VELOCITY FOR MERGING 

We follow B84 in estimating the critical velocity for merging of the tubes: if the 
magnetic energy per unit length of the reconnected flux exceeds the kinetic energy per 
unit length of the colliding tubes, the tubes coalesce. 

The magnetic energy Wm can be written in terms of an effective length Ai? e jj defined 

by 

AT? - 

It is important for the validity of our theory that -5 — e — be small, because we have solved 

for the structure of the reconnection layer only in an approximate way. It can be shown 

1 

that this ratio is just J(yo)y 4 , which is always less than 1 (see Figure 3). That is, only a 
small fraction of the fieldlines in the deformation layer are reconnected. 
The magnetic energy of the reconnected flux is 

W M = %,RAR eff = M+W** . (4.1) 

8tt jj 4 

The kinetic energy K per unit length of each tube is 

K = Meffvl = M* AC R*Bl(R) ^ 
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Using eqns (3.12), (4.1), and (4.2) we have 
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Wm 2 4 II 

-jT = -^R^I(yo), (4.3) 

with the criterion for coalescence > 1. 

A 

The function y^ 2 I(i)q) vanishes at both ends of its range (0,yo ;ma x), and attains a 
maximum value of about 0.35 at yo ~ 0.13. As R mg is expected to be large, Wm/K 
certainly exceeds unity at its maximum, so coalescence occurs over an interval in the 
collision parameter yo with endpoints given approximately by the roots of the equation 

We solve eqn. (4.4) by approximating /(yo) near the ends of its range using eqns (3.15) 
and (3.16). Near yo = 0, /(yo) ~ I m ax, and we find that coalescence occurs for 

^ > T2 9 7/2 R ' ( 4 - 5 ) 

or 

9I/2/-2 ,, 

M A c < ( 3/5 p ) 5/24 - (4-6) 

& ' taring 

The smallness of the lower limit of yo justifies the approximation / % I m ax- It can be 
shown froms eqns. (3.10) and (3.15) that the larger root of eqn. (4.4) is nearly j/o,mai 
itself, corresponding to the minimum velocity for reconnection, v cr a, given in eqn. (3.10). 
Putting these results together and using eqn. (2.10) we see that coalescence is expected 
for collision Mach numbers in the range M.AC,min to M.AC,max^ where 

M AC , ml n = — — ^TTPT, 4.7 

5 5 / 12 2 1 /8i?^ 12 



and 



M A C,max = ( ^^"V 24 = l.Q9M A C, ml n(^f /24 . (4.8) 



a 3 / 5 R mg a/i 



We will see in the next section that the range of Alfven Mach numbers for which coalescence 
occurs is between one and two orders of magnitude. 
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5. DISCUSSION 



In this section we critically review the main assumptions underlying our calculation. 
We then evaluate some of the formulae for parameters appropriate to the solar convection 
zone and discuss flux tube merging in that environment. Finally, we summarize the main 
results of the paper. 

5.1. Validity of Analysis 

We assumed in this paper that magnetic fluxtubes are thin and straight and have 
twisted field lines. The first two properties appear to be true of the tubes directly observed 
in the solar photosphere. Severely kinked or bent tubes are subject to strong magnetic 
tension restoring forces, so the assumption of straightness is probably a reasonably good 
one. There is no direct observational evidence that the field lines twist. However, magnetic 
fields which are concentrated by three dimensional flows possessing vorticity, as is the case 
in stellar convection zones, might well retain twist in addition to an axial component. 
Tubes which originate in the overshoot region at the base of the convection zone may also 
possess twist (Cattaneo, Chiueh, & Hughes 1990). The assumption of twist is central to 
the collision model. 

Following B84 , we considered direct collisions between antiparallel tubes. This is 
a maximally symmetric situation, and, as such, a singular one. Although extending our 
results to more general flux tube orientations is beyond the scope of this paper, the recon- 
nection process might not be too different if the tubes are not perfectly aligned, in which 
case some combination of the axial and azimuthal components would reconnect. 

The collision model assumes incompressible flows and slow collisions. If the relative 
velocities are comparable to the convective velocities deep in a stellar interior, and if the 
tubes are gas pressure dominated, the collisions are in fact likely to be slow. If the tubes 
are well coupled by friction to the surrounding gas, then the typical relative velocities of 
colliding tubes will in fact be much less than the mean convective eddy velocity. On the 
other hand, rapid relative motion is possible even with frictional coupling if the buoyancy 
and/or tension forces are strong enough (e.g. Moreno- Insertis 1986). Our conclusion that 
the amount of flux reconnected is a decreasing function of tube speed may change if we relax 
the assumptions that the collisions are slow, that the tubes are not greatly deformed, and 
that there are no important dissipation mechanisms other than resistivity. For instance, if 
a collision is strong enough to drive shocks in the tubes, there would be substantial viscous 
dissipation and the tubes might stay together long enough for extensive reconnection to 
occur. Likewise, tubes could be disrupted by violent, oblique collisions. However, for slow 
collisions, our flux pileup model should be generalizable to kinematic models other than 
B84 's. 

Finally, our description of the reconnection region is approximate, and based on the 
Parker- Sweet theory. We ignored additional field annihilation resulting from Ohmic diffu- 
sion. However, the contribution of the latter is expected to be very small. 

5.2. Application to the Solar Convection Zone 

In this section we apply our theory to the solar convection zone. Except where noted 
otherwise, our numerical values are derived from the convection zone model of Spruit 
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(1977). 



The global magnetic Reynolds number R mg scales as RB l . The denominator of this 



VP 



quantity is roughly constant throughout the convection zone and approximately equal to 
10 4 . Measuring R in units of 10 2 km and in units of 10 2 G we find 

R mg « 10 5 R 2 B 2 . (5.1) 

Note that the axial fieldstrength in the tubes is likely to be higher than B^; therefore R mg 
computed using the entire magnetic fieldstrength is significantly larger. 

We next compute the lower and upper limits v c i and v cu to v c for which our model 
predicts merging. The lower limit is given by eqn. (3.10) and using eqn. (5.1) it is given 
by 

vd ~ •34(^||)T2(p)-*cm5- 1 . (5.2) 
pii 2 

The upper limit v cu follows from eqns. (2.6) and (4.4), and is given by 

6.4, Bf . , , 

u » (5 - 3) 

The ratio of these two speeds is given by 

« 18.8( — - — ) 24 . (5.4) 

v c i fJ 2 a 

We see that the range of speeds for which our theory predicts coalescence is rather small; 
generally between one and two orders of magnitude. 

We now estimate v c i and v cu for parameters appropriate to the solar convection zone, 
and compare them with the characteristic fluid velocities v CO nv expected from mixing length 
theory. A range of velocities on different size scales is undoubtedly present in the convection 
zone, and as stated in the preceding section, the degree to which the tubes are coupled to 
the fluid motion is quite uncertain. Therefore, the mixing length velocity is only a rough 
benchmark. 

Halfway to the bottom of the convection zone the density is p % 6.4 X 10 -2 gm cm~ 3 . 
This leads to values of v c i and v cu of 1.4 and 26 cm/s, respectively, with all parameters in 
eqns. (5.2) and (5.3) having their fiducial values. On the other hand, v CO nv ~ 9.0 X 10 3 cm/,s. 
The magnetic fieldstrength B eq at which the kinetic and magnetic energies are equal at 
this depth is % 8000G; if we assume arbitrarily that ~ 0.1B eq} v c i is boosted to 4.8 
cm/s and v cu to 137 cm/s. We interpret the lower limit as essentially zero, while the 
upper limit is still a small fraction, no more than 1.5% in this example, of the convective 
velocity. 

At the base of the convection zone p % 0. 23 gm cm~ 3 , so v c i and v cw are reduced to .67 
and 12 cm/s, respectively. At this depth, v CO nv is only 86 cm/s. Increasing the magnetic 
fieldstrength to 0.1B eq makes little difference in this B eq is only 1460G; v cu would 

then be 16 cm/s. Small as this value is it is about 18% of the convective velocity, so a fairly 
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large fraction of flux tube collisions could result in coalescence by the processes included 
in our theory. 

Near the top of the convection zone, at a depth of 2 X 10 9 cm, p % 3.47 X 10~ 4 gm 
cm~ 3 . Then v c i and v cw are about 6 and 110 cm/s; again small compared with the 
estimated v conv of 2.1 X \0^cm/s. At this depth B e q ~ 4400(j, so assuming B2 ~ 4.4 
increases v c i to 14 cm s _1 and v cw to 360 cm/s, less than 1% of v conv . 

In summary, if the relative velocities between flux tubes are on the order of the 
convective velocity as obtained from mixing length theory, only a small fraction of collisions 
are slow enough to allow coalescence by reconnection. This fraction appears to increase 
monotonically with depth in the convection zone. 

Notice that the implied magnetic fluxes in our fiducial tubes are rather small; if the 
axial component of B is about a kilogauss, the flux is approximately that contained in 
a small photospheric flux tube. If the radii of flux tubes in the convection zone exceed 
100 km, eqns. (5.2) and (5.3) show that v r j and v cu are both reduced, although their 
dependences on R are rather weak. Stated in another way, if the tubes carry more axial 
flux (at constant axial fieldstrength) then v r j and v cu are both reduced, although their 
ratio is increased by Rmg 4 ■ 



5.3 Summary 

In this paper we have considered magnetic reconnection during slow collisions between 
twisted magnetic flux tubes. Our underlying assumption is that the tubes behave as 
described by Parker (1982), Moreno- Insertis (1986), and others: they are relatively long- 
lived, coherent structures and their motion is determined by frictional coupling to the 
background fluid, magnetic tension, and buoyancy. The frictional coupling may be strong 
but our model requires that it not be perfect, otherwise the dynamics of tube collisions 
would be determined entirely by the ambient flow and not at all by internal forces resulting 
from deformation of the tubes. Thus, for example, two tubes might intersect at an oblique 
angle and be dragged across each other: if magnetic tension is unable to snap them back 
and straighten them out, some forced reconnection is probably inevitable. This may occur 
in some situations. However, the analysis in this paper is restricted to flux tubes which 
are kept fairly straight (because of the flow, because of magnetic tension forces, or both) 
and which move semi-independently of the surrounding fluid. The tubes may be brought 
together because they are advected with the flow, through self-generated hydrodynamical 
forces (Parker 1979b), or some combination of them. Once the tubes are very close, the 
background flow must be essentially constant throughout the interaction region, so the 
dynamics of the collision should be governed entirely by internal forces. Although the 
overpressure associated with the deformation of colliding tubes tends to separate them, a 
sufficient amount of reconnected azimuthal magnetic flux will form a collar that keeps the 
tubes together. If the collision is too weak, there is so little deformation of the tubes that 
very little reconnection occurs, and therefore the tubes do not coalesce. If the collision is 
too strong it is so rapid that reconnection has very little time to occur, and again there is 
no coalescence. This is our basic picture, which we flesh out in the following paragraphs. 

We showed that Parker- Sweet reconnection of the magnetic field can occur at the 
interaction surface provided that the collision velocity exceeds a threshold which is ap- 
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proximately the geometric mean of the global Ohmic diffusion velocity in the tube and the 
Alfven crossing time of the tube, (measured with respect to the azimuthal field component 
only). The reconnection occurs through a very small amount of flux pileup at the edge 
of the interaction region, and only a small fraction of the flux in the interaction region is 
reconnected. There is a velocity which maximizes the amount of reconnected flux; it is set 
by balancing the driving, which increases with collision velocity, against the decrease in 
collision duration with increasing velocity, which reduces the time during which reconnec- 
tion can occur. If the collision velocity is less than the threshold, the fields are not driven 
hard enough to reconnect, but they will diffuse. Our results on reconnection are relatively 
insensitive to the values of the two parameters a and // (introduced in eqns. (2.6) and 
(2.8) respectively) which characterize the collision but are not well determined by the flux 
tubes model. We also argued in Section 5.1 that the essential features of this picture are 
preserved even if the tubes are not perfectly aligned. If the collisional velocity approaches 
the Alfven speed, our theory breaks down in many places, starting with the dynamics of 
the collision. 

We then discussed tube coalescence, arguing that merging occurs if the energy in 
the reconnected field exceeds the relative energy of the collision. This led to the result 
that merging occurs only for collision velocities within the range given by eqns. (5.2) and 
(5.3). We compared this velocity range with solar convection velocities estimated from 
mixing length theory, and found that they are always substantially smaller- although the 
gap narrows with depth in the convection zone. On this basis we argue that collisions 
between flux tubes do not generally lead to coalescence, although they might disrupt the 
flux tube edges. Again, our results are quite insensitive to a and //. The collision problem 
at high Mach number introduces effects we have not included, such as severe deformation 
and shock waves, and is probably amenable only to numerical analysis. 
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APPENDIX 

Our efforts to understand the restoring force in a collision of two parallel twisted 
magnetic flux tubes suggested that we investigate an analogous collision of two elastic 
solid cylinders. The usual heuristic comparison of Alfven waves traveling along a magnetic 
field line and mechanical waves on a string also leads to the comparison of MHD with the 
theory of elasticity. 

We assume an infinitely conducting magnetized fluid and consider the reaction of the 
system to a small perturbation away from an equilibrium state. The standard theory of 
elastic solids is likewise based on linearized equations describing the response of a system 
to a small perturbation from equilibrium. The equilibrium state is a minimum of the 
free energy, so the first nonvanishing term in the Taylor expansion of the free energy 
is generically of quadratic order. This leads to Hooke's law for an elastic solid and the 
usual linearized MHD equations for a magnetized fluid. We expect to need the equations of 
anisotropic elasticity, since the magnetic field is inherently anisotropic. Moreover, magnetic 
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fields possess directionality, a feature that is absent even in anisotropic elastic solids. This 
difference is expected to limit the ultimate utility of the analogy. 

The most general form of Hooke's law, applicable for anisotropic elastic solids (e.g., 
crystals), is 

Pki = Kklmdkdlim (-41) 

— * 

(Cf. Landau and Lifschitz 1970.) Here £ is the displacement of a volume element from its 
equilibrium position, which is assumed to be small, A is the elastic modulus tensor, p is the 
density, dj denotes partial differentiation with respect to spatial coordinate j, ' denotes 
differentiation with respect to time, and summation over repeated indices is implied. The 
elastic modulus tensor must satisfy the symmetries 

^iklm — ^ikml — ^kilm — ^Imik (-4-^) 

(Landau and Lifschitz 1970). 

We make a local approximation and assume a a uniform fluid at rest and a constant 
background magnetic field. The linearized magnetic induction equation can be integrated 
once in time to give the perturbed magnetic field SB in terms of the fluid displacement £ 
in coordinate notation as 

8Bj = B k d k Cj - Bjdkik . (A3) 

— * 

The linearized equation of motion for £ is 

Inp'ii = -di(B k 8B k ) + B k d k 8Bi = B k (d k SB z - d z SB k ) . (A4) 
We substitute SB as given by eqn. (A3) into eqn. (A4) to obtain 

4vrp li = -BjBkdidkij + B 2 didiii + BjBkdjd^i - BiBjdjd^i . (A5) 
We can rewrite this in a form resembling Hooke's law, 

47Tp£; = +B 2 S lk Sl m d k dl^ m -S lk BlB m d k dl^ m - BtBkSlrndkdl^rn + StrnBkBldkdl^rn . (A6) 

We then define 

:= BH ik 6 lm 

^^illm := —&ikBlB m — B z BkS{ m 

:= S im B h B x 

so that AW and A^ 2 ) are tensors with the symmetry properties of an elastic modulus tensor 
but the tensor A^ 3 -* lacks these properties. We note in passing that by orienting a coordinate 
axis along the uniform background magnetic field we reduce the number of nonvanishing 
components of the A^-* well below the maximum of 21 (cf. Landau and Lifschitz) allowed 
by the symmetry properties of the elastic modulus tensor. 

We now examine A^ 3 -* more closely. We can try to symmetrize it, defining 

\(3) ._ /x(3) w 4 

A iklm ■— \ A iklm A kilm A ikml A lmik)l^ 

= ^{SimBiBk + S km B\B l + 8uB m Bk + 8kiBiB m ) 
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We can now isolate the non-elastic part of the response, writing 

pi = + 2A< 2 > + AX^M^^U - BiBmdkdkU . (A9) 

If the last term on the right hand side is much smaller than the left hand side (i.e., small 
to the same order as other terms omitted in the derivation), then we have an approximate 
elastic modulus tensor describing the response of an MHD system to perturbations. In 
fact, we recognize that this term contains the scalar product of £ with B. One can show 
from the original equation of motion (A4) that this vanishes. Therefore, small amplitude 

— * 

motion is confined to planes perpendicular to B and is described by an elasticity tensor. 
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